Impurity scattering in the bulk of topological insulators 
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We study in this paper time-reversal 5-impurity scattering effects in the bulk of topological insu- 
lators (TI) in two and three dimensions. Specifically we consider how impurity scattering strength 
is affected by the bulk band structure of topological insulators. An interesting band inversion ef- 
fect associated with the change of the system from ordinary to topological insulator is pointed out. 
Experimental consequences of our findings are discussed. 
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A. Introduction 

Topological insulators (TI) are insulators characterized 
by a 7i2 topological number of the bulk band structure, 
and realize a new topological quantum state protected by 
time-reversal symmetry^. Theoretically, TI are insulat- 
ing in the bulk but possess gapless boundary states of he- 
lical Dirac fermions that give rise to interesting transport 
properties not realized in ordinary insulators (01). Novel 
excitations like Majorana zero modest have been pro- 
posed to exist in this class of materials and applications 
of the materials in spintronics have also been discussed 4 . 

Towards the applications of TI, a better understand- 
ing of the overall properties of the materials beyond their 
topological surface excitations is needed. For example, to 
utilize exotic boundary states of TI one needs to ensure a 
sufficient insulating bulk resistivity, which is a challenge 
at present due to the presence of bulk conductivity^— 
largely arising from the presence of impurities and defects 
in the materials^ 1 -. If the defects are external, which is 
the case for doped materials, one way to attain high bulk 
resistivity is to synthesize the materials with delicate bal- 
ance between donors and acceptor s] 12 ' 13 . On the other 
hand, we also want to understand how defects can affect 
the bulk conductivity through impurity induced in-gap 
bound states that lead to impurity band and effective 
narrowing of band gap. 

It was pointed out in Ref. |3| that point-like defects in 
TI do not give rise to topological protected zero-energy 
states. However, in-gap bound states can still be induced 
by (isolated) impurities*^— through conventional mech- 
anism. For example, depletion of wavefunction at the 
defect results in-gap bound states if one imagines that 
the impurity is formed by bending and shrinking the 
edge into a localized defect^. However, unlike the (real) 
edge states, these localized bound states are finite en- 
ergy modes because they suffer from symmetry allowed 
(self) interactions. In this paper we consider isolated 
time-reversal impurities described by 5-function poten- 
tials, and study how impurity scattering can depend on 
the special band structure of TPs. We shall show that 
impurity scattering can be enhanced in the bulk TI com- 
pared with 01 due to particular band structure associated 
with the materials. A criteria to search for materials with 
suppressed impurity scattering effect is given. An inter- 



esting sign-inversion effect in impurity scattering associ- 
ated with the band inversion in TI-OI transition is also 
pointed out. Experimental implications of these effects 
are discussed. 



B. Formulation 

We adopt the modified Dirac Hamiltonian in either two 
or three dimensions (2D/3D) as effective models for the 
topological insulator——, i,e, 



(1) 



k,.s 



where ip^Aipks) & re 2-component Dirac fields with mo- 
mentum and spin indices k and s =j, J., respectively and 



h 2 D(K s) = k x a y 



k y a x s z 



(2) 



for 2D topological insulators where cr's and s's are Pauli 



matrices, and rrik 



Bk 2 . a z = +!(—!) describes 



two contributing atomic orbitals to the topological in- 
sulator. For example, they represent s and p states for 
HgTe/CdTe system. The model can also be generalized 
to describe 3D topological insulators with a similar ef- 
fective Hamiltonian /i3£>(k, s) with = m — ^2 i Bikf 
where i = x,y,z and B.^s are parameters of the same 
sign. We shall for simplicity consider B x = B y = B ^ 
B z , corresponding to a common type of topological insu- 
lators in 3D with crystal structure R3m. 

For mBi < 0, the bands are ordered conventionally 
throughout the Brillouin zone and the system is a or- 
dinary (topologically trivial) insulator. The bands near 
k = are inverted due to strong spin-orbit coupling for 
mBi > and the system becomes a topological insula- 
tor. The energy eigenvalues of the bulk Hamiltonians are 



given by cjfe = ± y k 2 + (m — Bik 2 ) and are doubly 
degenerate (Kramers degeneracy). The time reversal op- 
erator is given by T = is v K with T 2 = —I, where K is 
the complex conjugate operator. 

In real materials there are other electronic bands and 
impurity scatterings exist in general between the Dirac 
bands and the other bands. To capture these effects 
we consider also inter-band scattering between the Dirac 
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bands with a model quadratic band described by 



ks v 



(3) 



where M and /i are chosen to have the same sign so that 
the quadratic band does not cross the Fermi level and the 
system remains a bulk insulator. M, fj, > (< 0) refers 
to a conduction (valence) band in this notation. 

We shall describe the impurity scattering by a single s- 
wave delta function potential sitting at the origin in this 
paper. In the four-band modified Dirac model we adopt 
here, the impurity scattering term within the topological 
bands can be written in orbital basis 1 — as 



Hji = iiS(x), 



(4) 



where u is a 4 x 4 matrix. For time reversal impuri- 
ties u should be Hermitian and time reversal symmetric, 
i.e. l~uT~ 1 = u. There are six bases for u, denoted as 



u 



(i=0,..,5). 



(I, a x , a\ a y s x , er a s a , a v s z ). 



To describe scattering between the topological and 
quadratic bands, we consider 



H 
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(*t( )V>„(0) + VU0)t*.(0)) , (5) 



([2]) is inadequate in describing the band-structure. Notice 
that the on-site Green's function matrix Go involved in 
Eq. (|6b[> is diagonal. This is a special feature of the delta- 
function impurity potential and would be absent if the 
impurity potential includes higher angular momentum 
components. 



C. g ± (u) around band-edge and enhanced impurity 
scattering 

Without loss of generality, we shall assume B > in 
the following. We note that in the limit of weak im- 
purity scattering, impurity-induced in-gap bound states 
can exist only if g ± (cj) diverges near the band edge. In 
usual (2D) semi-conductors the divergence in g ± (to) is 
originated from the band extremum associated with a 
quadratic band structure. In TPs, the band structure is 
more complicated and the band extremum may locate at 
a line of k values leading to stronger divergence behavior 
in g ± (oj). This phenomenon may arise in both 2D and 
3D TPs and will lead to enhanced impurity scattering 
compared with usual semiconductors. First we consider 
2D. 



where a = ±1 and s =t,|- ^s(x) is the Fourier trans- 
form of Ck s . We include only spin-independent scattering 
in writing down H12. 

The effect of impurity scattering is described by the 
T-matrix defined by 



(6a) 



T(u)=[I-UGo(u)) U, 



for the s-wave delta function potential, where U is the 
impurity scattering matrix and Gq(u>) — 57J k Go(k, uj) is 
the corresponding on-site matrix Green's function. In 
particular, the existence of impurity bound state is de- 
termined by the eigenvalue equatio n 19 i 20 



dot 



l-G (u)U 



0. 



(6b) 



For the modified Dirac model, the on- 
site matrix Green's function in the basis 
((1,1), (-1,1), (1,-1), (-1,-1)) (w.r.t.(a z ,^)) is 
given by 



.9 + M 
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where 



drijj / k 





D-l 



dk 



9 M - 



k 2 + 
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where Qd is the solid angle in dimension D = 2, 3. A is a 
high energy cutoff above which the effective Dirac model 



1. 2mB<l 

The integral ([5]) for g ± {oo) can be evaluated exactly 
in 2D. For 2mB < 1 the band minimum is located at 
k = with band gap \m\ as in usual semi-conductors. It 
is straightforward to show that in this case 



"MIm->m- 



-7T (lu ± to) In 1 



(9) 



where 5 ~ min(A, y |^|) is the momentum cutoff below 
which we may approximate to^ ~ to. We note that a 
logarithmic divergence appears in g ± (uj) near the band 
edge when \oj\ — > mT . Similar logarithmic divergence 
occurs for quadratic bands with Green's function given 
by 



g *(oj) = -2ttM\u 1 + 



A 2 



2M(/j - u) 



(10) 



The logarithmic divergence arises from the quadratic dis- 
persion of electron dispersion near the band edge and is 
responsible for appearance of impurity bound state for 
arbitrarily weak attractive scattering potential in 2D. 
Taking B > 0, we find that g + ^(uj) has logarithmic 
divergence at oj — > +(— )\m\ to — (+)oo for topological 
insulators (TI) (m > 0) and has logarithmic divergence 
at cj — y — (+)|to| to +(— )oo for ordinary insulators (OI) 
(to < 0), implying that the role of g + (~\oj) is inverted 
when the system changes from OI to TI. Microscopically, 
a system changes from OI to TI when there is an orbital 
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FIG. 1. Parameter range for band extrema ko- 



a ring of k values, k 
where ko = 



(k x ,k y ,0) with k\ 



k 2 



k 2 



mB — Region III is defined by 
B z > max(2^-,i?) with band extrema occurring at 
k = (0, 0, ±fco) where k 



l 



mB, 



The corresponding Green's functions are given by 



-47T 



k 2 dk 



(w ± m) 



ui ± m 



fc 2 



\frrv 2 



uj 2 tan 



const, as M 



Vfi 2 — <-<J 2 

(12a) 



inversion (2m = difference in energy between the two or- 
bitals forming the inverted bands) and the inversion in 
the role of g {uj) is a direct consequence of this orbital 
inversion. In other words, there is a direct correspon- 
dence between g^i(uj) 5oi( w ) with the same band gap 
\m\. 



2. 2mB>l 

In this case, the band extremum occurs at a line of 
momentum k with |fc| = ko- The corresponding bandgap 
is given by A = ^^"J^ 1 ■ Evaluating the integral, we 
find 



5 ± M|| W |- 



2B 



const. (11) 



We see that g (ui) has a stronger (inverse square-root) di- 
vergence at the gap edge compared with the case 2mB < 
1. This is a result of "dimension reduction" in the in- 
tegral when the band-minimum occurs at a ring of mo- 
mentum k with |fc| = ko. We notice that g + (ui) diverges 
to — (+)oo at ui — > +(— )A, with a weaker divergence 
(weighting factor) at uj — > — A. A similar situation oc- 
curs for g~~(uj) which diverges to +(— )oo at uj — > — (+)A, 
with a weaker divergence at uj — > A. The appearance of 
divergences at both uj = ±A is possible because of hy- 
bridization between the two atomic orbitals. We note 
that similar double-divergent behavior cannot occur for 
ordinary insulators which exist only at mB < 0. 



3D case 



We can study g ± {uj) for 3D TP s in a similar way. Us- 
ing the Hamiltonian h3n(k, s) with B x — B y = B ^ B z 
we find that the divergent behavior of g ± (u>) can be 
classified in three different regions as shown in Fig. [T] 
Region I is defined by 2m(B z , B) < 1 with band ex- 
tremum occurring at k = 0. Region II is defined by 
B > max (t^, _B Z ) with band extremum occurring at 
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0-2 

u}± Jb 
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dk. 
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k (uj± (m-Bkl)) 



k 2 + k 2 +m 2 k 



UJ' 
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(12b) 



A" 



-2vr 



E 



dk 

const, as M — s- A 



T fc °+Tr k(uj±(r 
dkz k 2 +kl 



B z k 2 )) 



(12c) 



where A = tttWvV AmB — 1, where B = B(B Z ) at region 

2 1 ts I 

II(III). Here S,5 Z are momentum cutoff as in Eq. ([9]). 
We notice that logarithmic divergence in g ± (ui — > ±A) is 
found for TI in region II which is absent in usual semicon- 
ductors which exist only at mB < 0. The sign- inversion 
effect in gf is also clear although there is no divergence 
compared with its 2D counterpart. 



4- enhanced impurity scattering 

The divergence behavior of g (uj) has strong implica- 
tions on impurity scattering behavior described by the 
T-matrix (|6a|) where we observe that impurity scattering 
can be strongly enhanced if g ± (ui) diverges. In particular 
TPs with non-diverging g ± (ui) , s are more robust towards 
impurity scattering and are better candidates for trans- 
port application. 

Thus (3D) materials belonging to region I and III are 
the preferred materials for transport application. For 
region II, due to the "dimension reduction" effect, diver- 
gence appears near the band edges at both ui — >• ±A and 
impurity scattering effects are much enhanced compared 
with usual semi-conductors. Notice that materials where 
all -B's are not equal do not suffer from "dimension re- 
duction" effect and are "preferred" also from the point 
of view of impurity scattering. 

We now apply our analysis to the three TI materials 
given in Ref. l2l[ where the low energy spectra around 
the inverted bands are given by 



Wfc — £fc ± 



B 2 k 2 



Ml 



(13) 
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material 


Bi2Se3 


Bi 2 Te 3 


Sb 2 Te 3 


band 


upper 


lower 


upper 


lower 


upper 


lower 


ko 


oA- 1 


0.09A" 1 


OA" 1 


0.04A" 1 


O.llA" 1 


0.08A" 1 


Region 


I 


II 


I 


II 


III 


III 



TABLE I. Band extrema ko and corresponding types in Fig. 
[T]of the inverted band for Bi2Se3 JL Bi2Te3 and Sb2Te3 (Band 
parameters are coming from Ref . I2H ) . 



where 
k 2 = kl 



dk 2 z + C 2 k z , M k = M + M x ki + M 2 k z and 
ky. Here we neglect the higher order term 
H3 which is relevant only for large k. Using the fitting 
parameters given in Ref. [2l] we classify both the upper 
and lower topological bands for the three materials in 
Table HI We see from TableUthat Bi 2 Se 3 and Bi 2 Te 3 both 
have type II band structure and Sb 2 Te 3 is the preferred 
material among the three as far as transport application 
is concerned. 



Impurity-induced in-gap bound states 

In this section we consider the formation of impurity 
induced bound states in the bulk where we shall illustrate 
the existence of an interesting band-inversion effect. First 
we focus on the intra-Dirac bands scattering. The scat- 
tering matrix respecting T can be written in orbital basis 
(the same basis as used in writing down Eq. ([7])) as 



u* 


U% 



-u* 




iti 

-U! 



Mil 
Ui 



Uj 

U22 



(14) 



where wn.22 £ K and ui,ui £ C. To see the qualitative 
effect of band structure and TI-OI inversion on bound 
state formation we consider the case of single scattering 
channels, i.e. when only one of the u's in Eq. (|T4"|) is 
non-zero. In this case it is straightforward to obtain the 
following eigenvalue equations 



(1- Ull g+(E)Y=0 
{I - U229- {E)) 2 = 
l-\u I \ 2 g+(E)g-(E)) 2 = 

{l-\u I \ 2 g+(E)g-(E)) 2 = 



(15a) 
(15b) 

(15c) 
(15d) 



for each non-zero it's where E is the bound state en- 
ergy. We shall analyse in detail the above equations for 
the case 2mB < 1 in 2D. The analysis can be easily gen- 
eralized to the case 2mB > 1 and to three dimensions as 
we shall see in the following. 

We start with inter-Dirac-band scattering. Since 
g + (E)g-(E) < (c.f. Eq. ©) in the gap region \E\ < 
|m|, the scattering matrix T has no divergence and there 



uu ^ (Eq. flgj)) 


wn > 


uii < 


m > 


No 


Yes 


m < 


Yes 


No 




"22 + (Eq. flgp)) 


"22 > 


"22 < 


m > 


Yes 


No 


m < 


No 


Yes 



TABLE II. Single channel intra-Dirac bands scattering in- 
duces in-gap bound state in the weak scattering limit, m is 
the "gap" parameter in the modified-Dirac Hamiltonian. 



is no bound state induced by inter-Dirac-band scatter- 
ing. The situation is very different for intra-band scat- 
tering where existence of in-gap bound state depends on 
the sign of Uug (E), and it is easy to see that a bound 
state which appears in the conduction/ valence band in 
01 (to < 0) state disappears in the corresponding TI 
(to > 0) state and vice versa, where we have again fixed 
B > for brevity. The results of bound state formation 
is summarized in Table [TTJ The doubly degeneracy of 
bound state solutions is a result of Kramer's degeneracy 
coming from time reversal symmetry. The existence of 
bound states in these cases is a natural result of an attrac- 
tive (repulsive) impurity in electron (hole) liquid, which 
is known to induce bound state in two dimensions for 
arbitrarily weak potential. The sign change in to just re- 
verses the conduction band to valence band or vice versa 
but the signs of u's are not reversed, resulting in the ap- 
pearance/disappearance of impurity bound states when 
to changes sign. 

We next consider impurity-induced inter-band scatter- 
ing between the Dirac-bands and an extra quadratic band 
described by Eq. ([5]). For simplicity we again consider 
single scattering channel with either u+(_ ) being non- 
zero. It is straightforward to obtain the eigenvalue equa- 
tions 

{l-vlg+(E)g*(E)) 2 = (16a) 
(l-v 2 _g-(E)g%(E)) 2 = (16b) 

where we again notice the Kramer's degeneracy of the 
solutions. Notice that unlike inter-Dirac-bands scatter- 
ing where g + (uj)g~ (u>) < 0, independent of the sign of 
to, g ± (E)gQ(E) depends now on the sign of m. The 
resulting bound state formation possibilities are summa- 
rized in Table ITTTl where we have considered the quadratic 
band to be either a conduction or valence band. Notice 
that the results are independent of the sign of the impu- 
rity potential since v a always enters the eigenvalue equa- 
tion as v 2 .. Again if we fix the scattering channel and 
the quadratic band, we find that the bound state ap- 
pears/disappears when to changes sign because the sign 
of g + (~> (E) changes upon orbital inversion. Our analysis 
can be extended easily to the case 2mB > 1. There is 
no band-inversion effect in this case since ordinary insu- 
lator exists only for mB < and impurity bound state 
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always exist because of the "double-divergence" behavior 
in g ± (uj). The situation for 3D TI's are similar. Band 
inversion effect exist only in Region I and weak impurity 
bound state always exist in Region II. 



v+ =/= (Eq. (US]}) 


M,/x > 


M,^ < 


m > 


Yes 


No 


m < 


No 


Yes 




v- ^ (Eq. QgbJ) 


M,|i > 


M, p < 


m > 


No 


Yes 


m < 


Yes 


No 



TABLE III. Induced in-gap bound states resulting from weak 
inter-Dirac-quadratic band scatterings. Quadratic bands are 
described by M,n > (conduction) or < (valence), and m 
is the "gap" parameter of the modified-Dirac equation. 



D. Summary 

We are now in the position to discuss and summarize 
our results. We show that impurity scattering effect in 
TI with band structure parameters in certain region are 
enhanced compared with usual semi-conductors. As a re- 
sult impurity bound states can form easily in these TPs, 
leading to enhanced bulk conductivity and reduced effec- 
tive band-gap. Impurity scattering effect can be reduced 
for TI's with band parameters in suitable region (region 
I and III in Fig. [T] for 3D TPs). Our analysis is applied to 
three materials Bi2Se3, Bi2Te3 and Sb2Te3 where we find 



that Sb2Te3 is the preferred material among the three as 
far as transport application is concerned. We note that 
a larger variety of TI materials exist nowadays and our 
analysis can be applied if their band structure is known. 

We also point out another interesting sign-inversion ef- 
fect associated with band inversion in TI-OI transition. 
The effect is most pronounced in 2D TI with 2mB < 1. 
Experimentally, the TI-OI transition can be controlled 
by gate voltage in InAs/GaSb quantum well system and 
the transport measurement can be carried out at various 
chemical potentials inside the band gap 22 i 23 . The band 
inversion effect can be observed through a change in dis- 
tribution of impurity bound state energies when the sys- 
tem changes from TI to 01 state. The effect is strongest 
if the impurities coupled preferably to one of the atomic 
orbitals. 

In conclusion, we study the effect of impurity scatter- 
ing in TPs in this paper. Our work is complimentary to 
previous work o that consider impurity bound state 
forming from continuous deformation of edge modes. We 
find that impurity-induced bound state formation de- 
pends strongly on the band structure of the TPs and 
impurity scattering can be suppressed in TPs with band 
structures in the correct region. Our result provides a 
guidance for TI material engineering which is useful for 
search of applicable TPs. 
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